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Abstract 

We investigate the possibihty that Us freezes as function of Nj within perturbation theory. 
We use two approaches - direct search for a zero in the effective-charge (ECH) (3 function, 
and the Banks-Zaks (BZ) expansion. We emphasize the fundamental difference between 
quantities with space-like vs. those with time-hke momentum. We show that within the 
ECH approach several space-hke quantities exhibit similar behavior. In general the 3-loop 
ECH (3 functions can lead to freezing for A'^^ ^ 5, but higher-order calculations are essential 
for a conclusive answer. The BZ expansion behaves differently for different observables. 
Assuming that the existence of a fixed point requires convergence of the BZ expansion 
for any observable, we can be pretty sure that there is no fixed point for Nf ^ 12. The 
consequences of the Crewther relation concerning perturbative freezing are analyzed. We 
also emphasize that time-hke quantities have a consistent infrared limit only when the 
corresponding space-like effective charge has one. We show that perturbative freezing can 
lead to an analyticity structure in the complex momentum-squared plane that is consistent 
with causality. 



1 Introduction 



The running of the strong couphng constant x{Q^) = agiQ"^)/!^ at large momentum 
transfers is well determined by the 1-loop perturbation theory result: 



and A is the QCD scale. Naively this leads to a divergent coupling at close 
to A (a "Landau pole"). However, higher-order terms, as well as non-perturbative 
effects are expected to alter the infrared behavior and remove the divergence from 
any physical quantity. 

It was suggested in the past, in various contexts [l]-[7], that infrared effects in 
QCD can be partially described using a coupling constant that remains finite at 
low-energies. A finite coupling in the infrared limit can be achieved in many ways 
without altering too much the ultraviolet behavior, which is well described by per- 
turbation theory. In the dispersive approach [4, 5], for instance, the "Landau-pole" 
is removed by power corrections, and therefore the infrared coupling is essentially 
non-pcrturbativc. On the other hand, it is possible [8] that the coupling constant 
freezes already within perturbation theory, due to a zero in the {3 function induced 
by the 2-loop or higher-order corrections. Such a perturbative infrared fixed point 
clearly occurs if the number of light fermions is just below Nf = 16| [8]. 

Phcnomenological studies, such as [1] and references therein, show that a coupling- 
constant that freezes at low energies may be useful to describe experimental data. 
But there is no general theoretical argument why QCD should lead to freezing. On 
the contrary, the general belief, which is supported by lattice simulations [9] and 
other approaches [10], is that for small Nf (and for pure Yang-Mills in particu- 
lar) there is no infrared fixed point. One appears only for large enough number 
of fermions. Thus there is some critical iV)"* (0 < Nf"^ < 16^) such that a fixed 
point exists only for Nf > Nj^'^^. The common lore is that below Nj"^ the theory is 
in a confining phase, with spontaneously broken chiral symmetry [11] and that the 
existence of an infrared fixed point is related to the restoration of chiral symmetry 



The question of the infrared behavior of the coupling constant immediately brings 
up the question of scheme dependence. The coupling constant is in general not a 
physical quantity, and only the first two coefficients of the (3 function are scheme- 
independent. It is clear, therefore, that the infrared stability of the couphng is 
scheme-dependent, and so is its infrared limit value, when it exists. Observable 
quantities can be used to define effective charges. These are useful, because contrary 
to arbitrary renormalization schemes, perturbative freezing of the coupling constant 
in physical schemes can provide some indication of the existence of a fixed point in 
the full theory. Clearly, the perturbative analysis can be trusted only if it leads to 
freezing at small enough coupling constant values. 




(1) 



where 




(2) 



[10]. 
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There are, in general, two techniques that have been used to study the infrared 
hmit of a generic observable within perturbative QCD. The first is a direct search 
for a zero in the f3 function, in a renormalization scheme that is 'optimized' to 
describe the particular observable [1, 12, 13, 14]. The second is the Banks-Zaks 
(BZ) approach [8, 2, 6, 15]. 

The 'optimized scheme' approach refers to a scheme in which one hopes that 
higher order corrections to the observable under consideration arc small. In this 
case the fact that one is using a truncated series for the observable and a truncated 
(3 function is hopefully insignificant. In particular, two schemes have been used to 
study freezing: the Principle of Minimal Sensitivity (PMS) [17], and the method of 
Effective Charges [16]. A detailed analysis of the infrared behavior based on the 
'optimized scheme' methodology was conducted in ref. [1] for total hadronic cross 
section in e+e~ annihilation {R^+^-), in ref. [12] for R^+e- and for the r lepton 
hadronic decay ratio {Rt)-, in ref. [13] for the Gross-Llewellyn Smith (GLS) sum 
rule for neutrino-proton scattering and in ref. [14] for the derivative of the hadronic 
decay width of the Higgs. 

The Banks-Zaks (BZ) approach [8] originates in the observation that in a model 
where the number of light quark flavors is just below N*^ = 16| (which is the critical 
value of Nf at which (3o changes sign), the perturbative (3 function (see eq. (3) below) 
is negative for very small values of the coupling constant, but due to the 2-loop term 
it immediately crosses zero and becomes positive. This perturbative infrared fixed 
point occurs at xpp — — 1/c = — /3o//?i, where (3i is the 2-loop coefficient of the (3 
fimction*. xpp — ^ as Nf approaches NJ. It was proposed [8, 15] to study higher- 
order effects on the fixed point for models with a varying Nf < Nf, by expressing 
xpp as a power series in the "distance" from Nf, i.e. in {Nf — Nf). It was later 
suggested by Stevenson [2] that this perturbative freezing of the couphng constant 
might be relevant for the real-world QCD with only two light fiavors. 

Recently the calculation of the 4-loop f3 function [25] has enabled Cavcny and 
Stevenson [6] to check the reliability of the BZ expansion for the location of the fixed 
point in the physical renormalization schemes defined through the effective charges 
of Re+e- a-nd of the Bjorken polarized sum rule, by calculating the O ((16^ - Nff^ 
order term in the corresponding BZ series. The authors of [6] have found that 
the relevant BZ coefficients are small for both observables they considered. They 
suggested that the small coefficients indicate that the BZ expansion indeed holds 
for Nf as low as 2, and therefore that the corresponding effective charges do indeed 
freeze due to perturbative effects in the real- world QCD. Caveny and Stevenson also 
investigated BZ expansion for the derivative of the Higgs hadronic decay width and 
for the anomalous dimension that is defined from the derivative of the /3 function 
at the fixed point. They found that the coefficients of these series are rather large. 
This may indicate that the BZ expansion breaks down at some larger Nf. This 
inconsistency was also one of the reasons we became interested in the subject. 

Another method for obtaining a finite coupling in the infrared limit from per- 

*The Nf dependence of /3i is such that it is negative for N^* < Nf < NJ, where Nf ^ 8.05 
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turbation theory, which attracted much interest recently, is the so-called "dispersive 
approach" or "analytic approach" [5, 4]. The idea is to construct a time-like dis- 
tribution by performing an analytic continuation of the running coupling constant, 
and then use a dispersive integral to define an effective space-like coupling that is 
free from spurious singularities such as a "Landau-pole" . 

The main purpose of this paper is to examine all the evidence for the existence of 
a perturbative fixed point in physical renormalization schemes and its dependence 
on Nf. We use both the 'optimized scheme' methodology and the BZ approach 
and study several different observables in order to get a global picture. We discuss 
relations between different physical quantities, such as the Crewther relation [21], 
and study their implication of perturbative freezing. We show and analyze the 
difference between quantities defined with space-like momentum and those with 
time-like momentum. 

This paper is composed of three main sections. The first (Sec. 2) is devoted to 
the 'optimized scheme' methodology and the second (Sec. 3) - to the BZ expansion. 
In third section (Sec. 4) we discuss infrared behavior of time-like quantities. The 
main conclusion are given in Sec. 5. 

The 'optimized scheme' part starts with a brief introduction on scheme depen- 
dence and optimized schemes in QCD (Sec. 2.1), followed by a computation and a 
discussion on the second renormalization group (RG) invariant (^2) for various time- 
like and space-like quantities (Sec. 2.2). In Sec. 2.3 we analyze the consequences of 
the Crewther relation regarding perturbative freezing. In Sec. 2.4 we examine the 
numerical proximity in p2 for several space-like quantities. In Sec. 2.5 we analyze 
stability of the fixed point in the Rg+e- effective charge in the 'optimized scheme' 
approach. The main conclusions of Sec. 2 are given in Sec. 2.6. 

The BZ section (Sec. 3) starts with a short introduction of the idea and the 
essential formulae (Sec. 3.1). We then present the BZ expansion in MS (Sec. 3.2), 
followed by an analysis of the BZ expansion for time-like quantities (Sec. 3.3) and 
a calculation of the BZ series for various quantities (Sec. 3.4). The consequences of 
the Crewther relation are given (Sec. 3.5) and finally we discuss the BZ expansion 
for the derivative of the P function at the fixed point (Sec. 3.6). The conclusions 
from combining the analysis of Sec. 2 and Sec. 3 are given in Sec. 3.7. 

In Sec. 4 we examine consistency of the perturbative approach and the "analytic" 
approach for describing the infrared region of time-like observables, such as Re+e- ■ In 
Sec. 4.1 we shortly review the dispersive relations between the vacuum polarization 
and Re+e-- In Sec. 4.2 we examine the analyticity structure of the D-function 
resulting from perturbation theory and show that it can be consistent with the one 
expected from causality only if the running coupling freezes. In Sec. 4.3 we examine 
the analytic perturbation theory approach. 
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2 Fixed Points from 'Optimized Schemes' 



2.1 Scheme Dependence and 'Optimized Schemes' 

We start by introducing the notation for the QCD /3 function, 

P{X) - = -pox'' - I3^X' - I32X^ = -l3oX^ (l + CX + C2X' + • • •) (3) 



where (3q is given in (2), the two loop coefficients is [22, 23]: 

/3i 1 



c = 



/So 4/3o 



102 - ^Nf 



(4) 



and higher-order coefficients C2, C3, • • • depend on the renormahzation scheme, and 
are given in MS by [24, 25]: 



C2 



16/?o 



2857 5033 ,, 325 ,,, 

Nf H iV? 

2 18 54 •'^ 



(5) 



/ 149753 



+ 3554(3) - 



1078361 6508 



162 



64/5o LV 6 

/50065 6472, \ , . 1093 
' H Cs Nf H TV? 

81 / ' 79Q / 



27 



+ 



V 162 



729 



(6) 



where („ is the Riemann zeta function ((^3 = 1.202, (5 = 1.03693). 

A generic physical quantity in QCD can be written in the form of an effective 
charge: 

x^^ = x (1 + rix + r2X^ + r2,x^ H ) (7) 

where x — ag/'K depends on the renormahzation scheme and scale. Specifying the 
expansion parameter x amounts to choosing all the coefficients of the fi function (cj, 
i > 2) and then setting the renormahzation scale. 

Consistency of the perturbative expansion (7), together with the RG equation 
(3), requires the invariance under RG of the following quantities for a given QCD 
observable [16]: 



p = ri - /3o In ( ^ 



(8) 



at first order, and 



P2 = C2 + r2 - Ti - ciri 

P3 = C3 + 2r3 + Arl + Cirl - 6rir2 - 2riC2 



(9) 



at second and third orders, respectively. Similar quantities can be defined at higher- 
orders [16]. 

At any finite order, a perturbative calculation has some residual renormahzation 
scheme dependence. In the infrared limit, the coupling constant in different schemes 
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can either diverge or freeze to a finite value. Therefore, the very existence of an 
infrared fixed point in a perturbative finite order calculation for any QCD quantity 
is scheme dependent. Of course, so is the value of the infrared coupling, when it is 
finite. 

Although theoretically any scheme is legitimate, in practice the choice of scheme 
is important even far from the infrared limit: there are schemes in which the pertur- 
bative scries (or the f3 function series) diverges badly and there are schemes in which 
a finite order calculation yields a good approximation to the physical value. In the 
following we shortly review two special schemes - the method of Effective Charges 
(ECH) [16] and the Principle of Minimal Sensitivity (PMS) [17]. These schemes 
are 'optimized', in the sense that the coefficients of the P function, as well as the 
renormalization scale, are set in a way suited to describe a specific QCD observable. 
As mentioned in the introduction, it was conjectured in the past [1, 12, 13, 14] that 
perturbative calculations in these schemes can be meaningful down to the infrared 
hmit, although, of course, the infrared values do not directly correspond to the 
measurable physical quantities, which are governed by non-perturbative effects. For 
Re+e- it was demonstrated [1] that if the experimental data are "smeared", they 
can be fitted by the perturbative result in the PMS^ scheme down to low energies. 

The ECH method [16] is based on using the actual observable effective charge 
as a coupling constant: xech = x^^. One way of achieving this is by choosing the 
renormalization scale and the renormalization scheme, i.e. the coefficients of the f3 
function, such that all the coefficients in eq. (7) are exactly zero*. It is easy to 
see that in this scheme, the coefficients of the j3 function are simply the invariants 
Pi listed in (9): cf*^^ = p2, c^^^ = ps, and so on. Thus, in order to find the value 
of the effective charge at the fixed point in a finite order calculation in the ECH 
scheme, one just looks for real positive solutions for the equation P^^^{x) — 0. For 
example, at the three-loop level one has: 

l + cx^^ + P2(x''^y ^0. (10) 

One further requirement is, of course, that x^^ at the fixed point will be small 
enough, so that the perturbative expansion will be trustworthy down to the infrared 
limit. It is clear that if c is negative and in particular if the c term dominates over 
the p2 term, there will be an infrared fixed point at approximately x^^ ~ —1/c — 
—Po/Pi- This is indeed the case in QCD with Nf just below Nf — 16| which is 
the critical value at which = 0. This is the starting point for the Banks-Zaks 
approach, to which we shall come back in Sec. 3. However, the physically relevant 
values for c are always positive (c is positive for Nf < 8.05). Then a real positive 
solution for eq. (10) can only originate from a negative p2- 

The PMS method [17] is based on using a renormahzation scale and scheme 
which is the least sensitive to a local change in the RG parameters. For instance, at 

^See Sec. 2.5 and Sec. 4 for a discussion on the applicability of the PMS/ECH methods to 

Re+e- ■ 

tThere are other choices possible at higher orders [12]. For instance at the three-loop order one 
can choose ri + = 0. These different possible schemes coincide at the fixed point. 
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the three-loop order, x^^ depends on two free parameters, which can be chosen to 
be the couphng x and C2. Thus the renormahzation scheme dependence of x^^ can 
be studied from the shape of the two dimensional surface of x^^ = C2). The 

PMS scheme corresponds to a saddle point on this surface at which 

dx 

and 



= (11) 



-t - ' '^^> 

The resulting equation, based on the condition P^^^{x) — 0, is [1]: 

| + cxpMS + 3|^P2- j^4ms = (13) 

and after solving the equation for .xpMs one substitutes it in (7) to get the value of 
x*^^ at the fixed point. Here, like in the ECH case, eq. (13) has a real and positive 
solution for a positive c and a negative p2- However, unlike the ECH case, here 
there is still a window for a small positive p2 (0 < p2 < c^/4), in which eq. (13) 
has a positive solution. In practice c is relatively small, and the conditions for the 
existence of a fixed point at the three-loop order are very similar in the ECH and 
PMS schemes. For both the ECH and PMS methods, the perturbative fixed point 
from the three-loop order analysis occurs at small coupling, provided p2 is large and 
negative. It was also found in [1, 12, 13] that the value of x^^ at the freezing point 
is somewhat lower in the PMS method than it is in the ECH case. 



2.2 The Second RG- invariant for Time- like and Space- like 
Quantities 

We saw that the existence of a perturbative fixed point in the ECH/PMS approach 
at the three- loop order depends crucially on the sign of p2- The fixed point found 
this way can only be considered reliable if it occurs at small enough value of the 
coupling - and this in turn depends on the magnitude p2- Therefore, the calculation 
of P2 is the first step in analyzing the perturbative infrared behavior of a QCD 
effective charge at the three loop level. 

In fig. 1 we present p2 bs a function of Nf for several QCD observables: 

1) The Bjorken sum rule for the deep inelastic scattering of polarized electrons on 
polarized nucleons, defined by 



dx 



gr{x,Q' 



1, , 

D 



1 



OiBj 



TT 



where 



aBj 

TT 



= XBj — x(l + kiX + k2X^ -\- 



(14) 



(15) 



with the perturbative result at NNLO from ref. [28]. 
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2) The Bjorken sum rule for deep-inelastic neutrino nucleon scattering, 

Fr{x,Q')-Fr{x,Q') 



/ dx 
Jo 



= 1 _ £^ 

~ 2 V TT 



where 



— = = + fiX + f2X^ H ) 

71 



(16) 
(17) 



with the perturbative result at NNLO from ref. [29]. 
3) The Gross-Llewellyn Smith sum rule (GLS) for neutrino proton scattering, 



/ dx 
Jo 



Fr{x,Q') + F,^^{x,Q' 



= 6 








TT . 



where 



OiGLS 



TT 



xgls = x{l + hx + kx'^ H ) 



(18) 
(19) 



with the relation to the coefficients of the polarized Bjorken sum rule given by 
[28]: 



h 



ho- 



rn 



'abc,]abc 



CfN, 



11 1 
144 ^ 6^^ 



Nf 



where the difference between the 3-loop coefficients of the two observables in 
(20) is due to the light-by-light type diagrams. 

4) The vacuum polarization D-function (not a directly measurable quantity) defined 
as the logarithmic derivative of the vector current correlation function n((5^), 
with a space-like momentum Q"^ = —q^ > 0, 

AttH I d'xe"^-^ {0\T{j^{x),f{0)} |0) = (g^ - q'gn^{Q') (21) 



DiQ') = Q 



where 



^,^_^,dUiQ^) 

OlD 




1 + ^ 
TT . 



TT 



= xd — x{l -\- dix + d2X^ -|- • • •) 



(22) 
(23) 



with the perturbative result at NNLO from ref. [30]. In (22) we ignored the 
contribution from the light-by-light type diagrams which is proportional to 
{J2fQfY- The corresponding diagrams contribute to x^) starting at three- 
loops: Axd = d2^x^ + 0(x'^) . We are interested in studying a purely QCD 
phenomenon, and therefore it is inconvenient to include these terms which 
involve also the electromagnetic interaction. However, it is still interesting to 
see whether this neglected contribution influences our conclusions concerning 
the infrared behavior. We will study this issue indirectly by comparing the 
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results for the Bjorken polarized sum rule to those for the GLS sum rule, since 
(cf. (20)) 

df^{h-k,)^^^. (24) 

5) The total hadronic cross section in e^e~ annihilation (again neglecting the light- 
by-light terms), defined by 



i?(.)^3(^gj 



i + — 

TT . 



(25) 



where 



— =XR^x{l + nx + r2X^ + • • •) (26) 

TT 

The perturbative coefficients of Re+e- can be related to those of the vacuum 
polarization D-function, by using the dispersion relation (see Sec. 4). The 
relations are [26]: 

ri = di 

rs = ^3 - 7r'/3o (^^1 + ^c) . (27) 

For our purpose, it is convenient to write the relations between the correspond- 
ing RG invariants pi defined in (9): 

= P?-l^'P'o (28) 

6) The T lepton hadronic decay ratio Rr, defined by [31], 

i^.^ ^(";^"-^"^^°"^y^ =3fl + ^V (29) 

where 

— = X-r — x{l + TiX + T2x'^ -|- • • •) (30) 
TT 

and where (7) indicates possible presence of photons in the final state. The 
perturbative coefficients of r,- are also related to those of the vacuum polar- 
ization D-function, as follows [31]: 

n = di-Poh (31) 

T2 = d2-{2di + c)Poh + Pol2 

n = ds- {U2 + 2dic + C2)(3oh + Udi + \c] (3^2 - Poh 



2 

where h = -19/12, I2 = 265/72 - and I3 = -3355/288 + 197rVl2 
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7) The static potential, defined by [32] 



V(,-) = -c'Jl^ (32) 

where is the three-momentum squared, corresponding to the spatial sepa- 
ration r between the quark and the anti-quark, and where 

— = xv — x{l -\- Vix + V2x'^ + •••). (33) 

TT 

The static potential was recently calculated [32] up to order 0{al). 
8) The derivative of F^, the Higgs hadronic decay width, defined by [33] 

an _ ld\n{TH/MH) ^ , , i, 2 , x /q^in 

— = xh^-- 2 ^x{l + hix + h2X +■■■) (34 

TT 2 a in Mh 

where Mh is the Higgs mass, which is assumed to be much larger than the 
quark masses, xh was recently calculated [33] up to order 0{a^). This means 
that p3, the 4-loop coefficient of the ECH (3 function is now available. 

The observations from fig. 1 are: 

a) There is a clear distinction between the time-like quantities and the space-like 

quantities. 

b) There is a surprising numerical proximity between for several different space- 

like quantities: this includes the vacuum polarization D-function, the GLS 
sum rule and the polarized and non-polarized Bjorken sum rules, but not the 
static potential. The proximity is particularly evident for Nf < 5. This issue 
and its relation to the ideas of ref. [18] are further discussed in Section 2.4. 

c) For the space-like quantities (except the static potential) p2 becomes positive 

for Nf ^ 4, and thus according to the ECH/PMS approach at this order there 
is no fixed point for these values of Nf. 

d) The static potential behaves differently from the other space-like quantities, 

becomes positive already for Nf ^9. 

e) At low Nf, P2, P2 and are numerically very different from one another. For 

larger Nf they become closer. At Nf > 8, the p2 values for the three effective 
charges are large and negative and very close to each other^, indicating that 
they freeze to similar values (sec fig. 2). 

f) For the time-like quantities p2 stays negative down to Nf = or 1, indicating a 

possible infrared fixed point according to the ECH/PMS approach. 



^This is due to their relation through an analytic continuation as discussed in Sec. 4. 
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g) P2 (p2 for the derivative of the Higgs decay width) behaves differently from the 

others: it is not a monotonically decreasing function of Nf, and it is negative 
and large for any Nf. 

h) The second coefficient of the (3 function in MS (02) is not close to those of the 

effective charge schemes (^2)- 

Let us briefly discuss the question of freezing of xh, the effective charge deflned 
from the derivative of the Higgs decay width, as we believe it can teach us a general 
lesson. Naively, the fact that is negative and large means that there is an infrared 
fixed point at a rather small Xh value. For instance, for Nf = 3, ~ —57.6, and 
thus the ECH [3 function (see eq. (10)) has a zero at x^^ ~ 0.15. Now, since the 
value of Xff^ is small, one could further conclude that the perturbative analysis is 
reliable. But is this really the case? This issue was discussed in detail in ref. [14], 
where is was conjectured that this fixed point is spurious, and that the perturbative 
series breaks down at NLO in this case. In general, the fact that p2x'^ becomes 
equal to the leading terms 1 + cx does not immediately imply breakdown of the 
perturbative series^. It is still possible that p^x^ (and maybe a few higher order 
terms) will be smaller than p2x'^, while the asymptotic nature of the series will take 
over at some higher order (for a recent review see [41]). Since the next order term 
for Xh is now available [33] this question can be clarified explicitly: it turns out 
that indeed the 4-loop ECH (3 function —(3qx'^ {l + cx + pf^a:^ + pfx^^ has a real 
and positive zero only for Nf ^ 14. For Nf ^ 13 the term p^x^ turns out to be 
larger than p^x^ before p2X^ becomes larger than the leading terms 1 + cx. This 
confirms that the conjecture in ref. [14] that the fixed point in xh (for a small Nf) 
is a spurious one and that the series breaks down at NLO. It is clear that all-order 
resummation is essential in this case. The general lesson is that one should exercise 
extreme caution when looking for a fixed point in a finite order calculation. The 
NNLO analysis of the existence of a fixed point in the ECH/PMS schemes is based 
on the assumption that \p3Xfp\ < \p2\, which may turn out to be wrong. 

Now we proceed to consider other quantities: in fig. 2 we present the ECH 
value of the effective charges at freezing obtained through solving eq. (10), for the 
D-function and for the related time-like quantities: Re+e- and i?^. Of course, this 
calculation can only be meaningful if p2 < 0. Prom (10) it is clear that x^p diverges 
as p2 — > 0", corresponding to A^/ ~ 4 for xd and to A^/ ~ 1 for Xr- 

The Xd result in fig. 2 actually represents very well also the results for the 
other space-like quantities (except the one for the static potential xy). The reason 
is transparent from figs. 1 and 3 - the latter showing the differences between the 
values of p2 for various quantities and pf : the differences between p^, pf-', p2^ and 
P2^^ are rather small in the region Nf < S, where the three loop contribution is 
important for freezing. For higher value of Nf, freezing is induced by the two loop f3 
function, which is invariant, and the differences in p2 almost do not alter the value 
of the effective charge at freezing. 

^If that was the case, then 'perturbative freezing' would have been a meaningless term. 
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In fig. 2 we show the PMS result for the D-function only. The PMS results 
for the effective charge at freezing are somewhat lower than the ECH result (in 
accordance with [1, 12]) but the general picture is the same. 

Purely perturbative effective charges, at any order in perturbation theory, can in 
principle diverge in the infrared, independent of whether or not the full theory has 
an infrared fixed point. A priori, it is also possible for different effective charges to 
have a totally different perturbative infrared behavior for a given Nf. In particular, 
it is possible that there will be a "Landau-pole" in one effective charge Xa, while 
another effective charge xi, will exhibit perturbative freezing. Still one could write a 
power expansion of x^ in term of Xa (and vice- versa) [18, 20]. One would then expect 
such an expansion to be divergent. However, in practice the numerical proximity 
between the p2 coefficients for several different space-like effective charges suggests 
that the expansion of Xa in terms of Xh is at least close to being convergent (See 
Section 2.4). It seems that the different space-like effective charges considered above 
(excluding xy) are so closely related that perturbative freezing could only occur for 
all of them together or - for none. 

2.3 Freezing and the Crewther Relation 

The only example where one can explicitly verify our conjecture about simultaneous 
perturbative freezing of different quantities is the Crewther relation [21]. This is an 
all-order relation between the vacuum polarization and the polarized Bjorken sum 
rule. In terms of effective charges it can be written as 



where the /3 function is defined in (3), S{x) is a power series in the coupling constant 



and Si depend on Nf and Nc- Writing the effective charges as power series (eqs. 
(15) and (23)) one obtains a relation between the coefficients of XBj and di of xd, 
as follows: 



Prom the knowledge of the 3-loops coefficient of the Bjorken sum rule [28] and the 
vacuum polarization [30] we can obtain Si and 5*2 : 



XBj - Xd + -CfXBjXD 



-P{x)S{x) 



(35) 



S{x) = Si + S2X + S^x^ + ■ ■ • 



(36) 





+ I2C3 



(38) 
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and 



221 629 38 397 f,nr r m M C^o^ 

where 5*2 is scheme dependent and is given here in MS . 

The term on the r.h.s. of eq. (35) is scheme-invariant, but P{x) and S{x) are 
separately scheme dependent. If xd has a perturbative fixed point x^^, then it is 
convenient to write the r.h.s. of (35) in terms of x^. I3[x^^) = and so the r.h.s. 
vanishes at xd = x^^. Therefore xbj also freezes perturbatively, leading to the 
original conformal Crewther relation: 

The argument works, of course, in both directions, i.e. if the Bjorken effective charge 
freezes, then the D-function will also freeze to the value: 

^FP 
1 — l^fX^j 

Note that (35) allows a situation in which both xd and xsj diverge in the infrared 
limit. 

We shall use the Crewther relations in the following section, where we study the 
numerical proximity of p2 between various space-like quantities, and also later, in 
conjunction with the BZ approach. 



2.4 Numerical proximity of p2 for different space- like quan- 
tities 

In this section we study the numerical proximity between the p2 invariants for the 
various space-like quantities (fig. 1) mentioned in Section 2.2. The values of pf'', 
P2, pf^ and P2^^ are given in Table 1, for iV^ = 3 and Nf < 7. In fig. 3 we show 
the difference between the values of p2 for the various space-like quantities and pf*. 
The vertical scale is enlarged here by a factor of 10 with respect to that in fig. 1. 



The numerical proximity between p2 for the GLS sum rule and p2 for the Bjorken 
polarized sum rule is obvious, as they differ only by the light-by-light diagrams at 

3-loops level, giving rise to a small" contribution, proportional to Nj. This difference 
of course vanishes exactly for Nf = 0. Therefore, we focus on the three remaining 
quantities: the vacuum polarization D-function, and the polarized and non-polarized 
Bjorken sum rule. 

II It is interesting that this difference becomes important in the framework of the BZ expansion, 
as we discuss in Sec. 3.4. 
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Nf 


Bj 











17.812 


17.924 


15.956 


17.924 


1 


13.557 


13.747 


12.344 


13.334 


2 


9.371 


9.605 


8.702 


8.779 


3 


5.237 


5.475 


5.007 


4.236 


4 


1.134 


1.330 


1.223 


-0.322 


5 


-2.9(39 


-2.869 


-2.691 


-1.935 


6 


-7.114 


-7.177 


-6.798 


-9.657 


7 


-11.361 


-11.669 


-11.186 


-14.562 



Table 1: p2 for various space- like quantities for Nf <7 



One's initial guess is to suspect that the similarity between p2 for the Bjorken po- 
larized sum rule and p2 for the vacuum polarization D-function is due the Crewther 
relation (see Sec. 2.3). We show below, however, that this numerical agreement is 
a second "miracle" , on top of the Crewther relation. 

The numerical agreement between pf^ and the rest of the space-like quantities 
at low Nf is not as good: Bit Nf — the difference is about 10%, vs. 0.6% relative 
difference between p^-' and p^. Nevertheless, for larger Nf pf^ is quite close to the 
others. 

We now investigate further the numerical proximity of p^'^ and p^, for two rea- 
sons: first, the numerical agreement in this case is remarkable for < A^/ < 7 (cf. 
Table 1)**, second, it is interesting to see to what extent this numerical proximity 
is related to the Crewther relation. 

Prom (9) and (37) we get: 

-P2= \CfPoSi + IPoS2 - Id.PoSi - ^P'oSl + lcCf (42) 

where di is defined in (23) and 5*1 and 5*2 are given in (38) and (39). Note that both 
di and 5*2 depend on the renormalization scheme and scale, but in such a way that 
pf-' — P2 is scheme and scale invariant. Substituting the three- loop expressions into 
(42) we obtain a complicated function of Nc and Nf, with no clue that pf-' — pf is 
small. To make things simple, we first consider the case Nf — 0: 



Bj D 
P2 - P2 



_ ( 1043 869 55 121 

iV/=o - ^4" 4752 " 108 ^=^"¥^^ + "9"^^^ ^^^^ 
23053 253, 55, 



Numerically, we have: 



Bj D 
P2 - P2 



19008 36 ^ 6 



~ 0.02966A^,^ - 0.15542 (44) 

Nf=G <^ ^ ' 



**The agreement is not as good for larger Nf values, for instance, for Nf = 14, p^'' = —80.2 
while P2 = —68.8, about 14% difference. This is to be compared with a difference of 1% to 3% for 
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The "miracle" is in the numerics! In (43) one finds all the irrational numbers that 
enter the three loop calculation (^3, ^| and (^5. Each of the terms separately is of 
order 1 to 10, but they combine to give a tiny sum. 

It is important to note that there is nothing special in the Nf = case considered 
above. In order to get a more general view of the relative magnitude of difference 
P2'' — P2 : we consider the normalized difference TZ defined by: 



7^ = 



1 Bj 
W 




|P2 1 





(45) 



IZ is plotted in fig. 4 as a function of Nc for various values of A^/, and in fig. 5, 
as a function of Nf, for various values of Nc- In general, TZ is of order 1%! The 
only occasions where the relative difference TZ is not small (the peaks raising above 
TZ = 0.1 in fig. 4 and 5) is when both pf"' and P2 are close to zero. Then they may 
even have opposite signs, leading to 7^ = 1. 

We conclude that the numerical proximity between p2 and p^-' is not a direct 
consequence of the factorization implied by the Crewther relation (r.h.s. in eq. (35)), 
but of the particular numerical coefficients. While the numerical proximity of ^"^ 
and P2'' is well understood, we do not know of any reason why p2^ is close to p^"'. 
It is tempting to think that there is a deeper reason for the numerical proximity, 
and that higher order ECH coefficients (pi for i > 2) for different quantities are also 
close. In this respect it will also be interesting to know the fundamental reason why 
P2 is so different than p2 for the other space-like quantities. 

Next we discuss the relation between the assumption that pi for different observ- 
ables are close to one another, and the work by Brodsky and Lu [18] on commensurate- 
scale relations between observables. In [18] (see also [20]) it was suggested to express 
one effective charge (xh) in terms of another (xa), and then choose the scale of Xa 
according to the BLM criterion [19]. In [18] it was found that the coefficients in 
the expansion Xh = Xh{xa) are both much simpler than the ones in some arbitrary 
scheme (like MS), and are numerically small. Let us see what are the conditions 
for the coefficients in such an expansion to be small. We start with the expressions 
for two generic effective charges in some scheme: 

x^ = x + rlx^ + rgX^ H (46) 

xi, = X + r^x^ + r2X^ + ■ ■ ■ 

and express Xh in terms of Xa- 

Xb^Xa + mixj + m2Xa H (47) 

where depend on rf and r\, for instance: mi — r\ — r^, m2 — r\ — r2 — 2r1{r1 — r\) 
Note that are, by definition, invariant with respect to the choice of the inter- 
mediate renormalization scheme. The next step is to use the definitions of the RG 
invariants (9), and express the coefficients of (47) in terms of pi and p\. mi is just 
the difference between the values of the first invariant p (8) for the two effective 
charges: 

mi = - = p*" - p" = A, (48) 
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and higher-order coefficients can be expressed entirely in terms of A and of the 
higher-order invariants and p^, and depend mainly on A and on the differences 
of the coefficients pi. We define Si — — , and then, 

1712 = 52 + A.'^ + cA (49) 
ma = ^53 + A(p^ + 252) + ^cA2 + A^ 

A can be tuned by changing the scale of Xa, while the scale of xi, is kept fixed. In 
particular, there are choices of scale for which A is small, such as the leading-order 
BLM scale [19, 18], which eliminates all f3o terms from A, or, simply the choice 
A = A small A is, however, not enough to guarantee small higher-order rrii 
coefficients. The latter will be small only if 5i are also small. 

One practical conclusion from this discussion is that relating xd, xbj, xqls and 
xpi to one another as suggested in [18] will probably lead to more accurate results 
than some generic scheme. However, relating any of the above to Xy is disfavored. 

This concludes our discussion on the numerical proximity of the p2 coefficients 
for space-like effective charges. Next, we briefiy discuss the possibility of applying 
the 'optimized-scheme' approach to study the infrared limit of time-like effective 
charges. 

2.5 The reliability of the "fixed point" in Re+e- from ECH 

In this section we study further the 'optimized scheme' approach applied directly to 

the Re+e- effective charge, along the lines of rcf. [1, 12]. We consider specifically the 
case Nf = 3 at the three- and four-loops order, and try to estimate the reliability 
of the ECH analysis in the infrared. A deeper study of time-like quantities in the 
context of freezing is postponed to Sec. 4. 

The 'optimized scheme' approach was found to be unreliable when applied di- 
rectly to time-like quantities in another context, by Kataev and Starshenko [34]. 
They found that the ECH and PMS methods for estimating the next term in a 
series* when applied directly to two- and three-loop series for time-like quantities 
such as Re+e- , do not predict the correct structure of the terms that result from the 
analytic continuation (the situation is similar for Pade approximants) . The solution 
of [34] is to predict the coefficients of the space-like D-function and use the exact 
relations (27) to obtain the coefficients of the time-like quantity. One will naturally 
expect that if PMS and ECH fail in predicting the next terms when applied directly 
to the time-like quantities, they should not trusted for studying the infrared limit 
of these quantities. Nevertheless, it may still be instructive to see what one obtains 
in this approach. 

^^An advantage of the first over the latter is that in the first the scale does not depend on the 
number of light flavors, and thus observables cross the quark thresholds together. This issue is 
discussed in detail in ref. [18]. 

*These methods arc based on the assumption that the next, uncalculated coefiicient in the 
ECH/PMS /3 function, is close to zero 
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Considering eq. (28) one finds that the tt^ terms, that make the coefficients of 
the ECH /9 function of (pf ) different from those of xu {pf)-, are numerically 
significant and negative. This is the reason why the analysis in Sec. 2.2 indicates 
an infrared fixed point for the xr and not for xd for Nf < 4. 

Taking as an example QCD with Nf = 3, we examine the ECH /3 function for 
the D-function and Re+e-' 

p^-pox'^{l + cx + p2X^ + P3X^ + ■■■). (50) 

The two loop coefficient is invariant: c = 1.778. The three loop coefficients are 
P2 c::^ 5.23 and ~ —11.42, and the four loop coefficients are 

pf ~ -33.39 + 2(^3 (51) 

and 

pf ~ -181.98 + 24 (52) 

where ^3, the MS four loop coefficient of the D-function, is unknown. 

Next, we consider different approaches to predict 0^3, and thus and p§. In 
the ECH/PMS approaches for predicting the next term in a perturbative series 
[16, 17, 34] one assumes pf ~ in order to obtain a prediction for ^3. Therefore, 
one cannot use an ECH/PMS prediction to calculate pf . Fade Approximants (PA) 
[36] applied in MS predict ^3 = 24.75 (using the [1/1] PA) and ^3 = 16.49 (using 
the [0/2] PA). We note that these predictions are close to one another, and are also 
consistent with the ECH/PMS assumption pf ~ 0, which leads to ^3 = 16. 7^ One 
could alternatively try to use either PMS/ECH directly for the Rg+e- (i-e. assume 
P2 = 0, instead of assuming pf = 0) or apply the PA's method directly for the 
Re+e~ series. However, as we already mentioned, the resulting predictions do not 
agree between the different methods and do not contain the correct tt^ terms [34] 
and are therefore not rehable. 

We conclude that pf ~ and pf ~ —150. Thus, while the question of whether 
the D-function effective charge freezes for Nf = 3 remains open, it is quite clear 
that the naive 'optimized scheme' approach predicts that Re+e- does: pf is large 
and negative. Nevertheless, if the value of xr at the fixed point is calculated (as the 
zero of the four-loop ECH /3 function) using a reasonable guess for ^3, one obtains 
a x^^ ~ 0.19 ± 0.03 * which is compared in fig. 2 with the three-loop ECH result 
x^^ ~ 0.38 (or the three-loop PMS results x^^ ~ 0.30). We therefore conclude that 
the ECH (or PMS) methods fail in predicting the infrared limit of the perturbative 
result in this case. The reasons for this will become clear in Sec. 4. 

A remark is in order concerning the Higgs decay width effective charge xh, which 
is also a time-like observable, and therefore contains tt^ terms that result from the 

'^The authors of [34] obtained a sightly larger value: ds = 27.5. The difference is due to the 
additional assumption taken in [34] that C3 = - the 4-loop coefficient of the /? function was not 
known then. 

■I-The only real solution of the equation /J^^h _ q a,t 4-loops (with Nf = 3) is pretty stable: it 
changes in the range 0.16 < x^^ < 0.23 for -20 < ^3 < 60. 
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analytical continuation. In principle, the problem discussed above should appear in 
this case as well. However, contrary to the case oixn or the numerical significance 
of these terms in the xh series is rather small (compared to other contributions at 
the same order) and therefore our previous conclusions concerning xh hold. 

2.6 Conclusions 

We analyzed here the freezing of the QCD effective charge for various quantities by 
looking for a zero in the corresponding ECH (or PMS) (3 function. 

We found a clear distinction between time-like and space-like quantities. For 
several space-like quantities we found that the numerical values of the second RG 
invariant are quite close, especially for a small Nf. This suggests that the different 
quantities are closely related. It is tempting to conjecture that this will be reffcctcd 
in numerical proximity of pa and higher-order coefficients of the corresponding ECH 
(3 functions. We also expect that when perturbative freezing occurs, it will occur 
together for the various quantities. This is clearly true for the D-function and the 
Bjorken polarized sum rule, due to the Crewther relation. For other quantities, such 
as the Bjorken non-polarized sum rule, this is just a conjecture. 

For the space-like quantities (except for the static potential) it seems possible 
that there is a fixed point for Nf ^ 5, while there is no indication of freezing below 
Nf ^ 5. Absence of a perturbative fixed point at the three- loop level does not 
necessarily mean that it does not exist. It simply means that one should consider 
higher-order correction to answer this question. On the other hand, presence of a 
fixed point at the three-loop level, does not guarantee that it will persist at higher 
orders. 

In Sec. 3 we look for more clues for or against the relevance of perturbative 

freezing, by studying the Banks-Zaks expansion of these quantities. 

For time-like quantities, a naive application of the PMS/ECH approach indicates 
a perturbative fixed point even when the corresponding space-like quantity does not 
freeze. On the other hand we identified an instability of the predicted effective 
charge value at freezing. This result is a first indication of the inconsistency of the 
approach, as discussed further in Sec. 4. 

3 The Banks-Zaks Expansion Approach 
3.1 The BZ fixed point 

We start this section by summarizing the basics of the BZ expansion [8, 6, 15] for 
the location of the fixed point. 

As is mentioned in the introduction, the BZ expansion is an expansion in [Nf — Nf), 
i.e. the "distance" from the critical value Nf — 16^ down to lower values of Nf. It 
is convenient to use the expansion parameter [15, 6]: 
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Since oo is linear in Nf, it is straightforward to rewrite the (3 function coefficients 
Pi — PoCi as polynomials in qq, obtaining (cf. (2), (4)-(6)): 

c = + ci,o = + — 54 

and 

1 , , 

C2 = C2 -1 h C2,0 + C2,iao (55) 

ao 

C3 = C3 _i h C3,o + C3,iao + C3,2ao (56) 

ao 

where explicit expressions for Qj- in MS are given in ref. [6]. The next step is to 
solve the equation (3{x) = 0, yielding an expression for xpp. With the 2- loop (3 
function, one obtains: 

Prom (57) it is clear that xpp is asymptotically proportional to oq, and therefore we 
look for higher-order solutions xpp, in the form of a power expansion in ao'- 

xfp = ao{l + Uiao + U2al-\ ). (58) 

Assuming that a fixed point exists, one then substitutes (58) in the equation j3{x) = 
0, finding that Ui is fully determined by the 3-loop coefficient of the f] function, U2 
- by the 4-loop coefficient, etc. The formulae for the first three Wj-s are: 

Ul = Ci,o + C2 -1 

U2 = cI q + 2cl _i + 3ci,oC2 -1 + C2,o + C3 _i 

^^3 = 3c2,oCi,o + 4C2,0C2 -1 + 4Ci,oC3 _i + C^ q + 6cI qC2-1 

+ 10Ci,oC2 _i + C4 -1 + 5C2 _lC3 _i + 5C2 _i + C2,l + C3,o (59) 

Ul, U2, ■ ■ ■ depend on the renormalization scheme. One may be interested in studying 
the fixed point in some other scheme, related to the original one by 

y = x{l + rix + r2X^ + •••). (60) 

In particular, one may be interested in the fixed point in some physical scheme, in 
which case y is an effective charge corresponding to some measurable quantity. In 
order to obtain the appropriate expansion for the location of the fixed point ypp in 
the new scheme, one first writes the coefficients as polynomials in gq: 

n = Vifi + Vi^iao H h Tibial- (61) 

Next, one substitutes (61) and (58) in eq. (60), obtaining 

ypp^ ao{l + Wiao + W2al-\ ) (62) 
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where 



wi = ui + ri,o 

W2 = U2 + 2rifiUi + r2,o + n^i 

W3 = «3 + 2ri,oM2 + 2ri,iMi + 3r2,o«i + ri,o«i + ^^2,1 + ra.o (63) 

It is important to note that the coefficients Wi for a given effective charge y are free 
from any renormahzation scheme ambiguities [15], as scheme dependence cancels 
out between the Ui and terms. This is expected on general grounds, since both 
Ufp and oq in eq. (62) are physical quantities. This invariance can also be under- 
stood considering another possibihty of obtaining the BZ coefficients: as a first step 
one calculates the pi renormahzation scheme invariants (9) of the required effective 
charge. One then writes the coefficients pi as series in Oq, similarly to (55) and (56), 
with the only difference that here also pi^i terms appear, for instance: 

P2 = P2,-l — + P2,Q + p2,iao + ^2,2^0 (64) 

ao 

Finally, one calculates the BZ expansion for the value of x^^, directly from the ECH 

(3 function f3{x^^) = —Pq (x^^^ 1 + cx^^ + p2 (x^^^ + p^ (x^^Y + • • • , similarly 

to the way eq. (59) was obtained for the MS (3 function. The resulting BZ series 
is exactly equal to the one obtained in (63) using the BZ expansion in MS as an 
intermediate step. 



3.2 The BZ Expansion in MS 

Before using the BZ approach to study perturbative freezing for effective charges we 
consider the expansion in MS , where the coefficients of the f3 function are given by 
eqs. (2) and (4) through (6). We are interested the range < Nf < 16|, in which 
there is asymptotic freedom. We note that c is negative (leading to a fixed point at 
the two loop order) only for Nf > 8.05; C2 is negative only for Nf > 5.84, and C3 is 
never negative. This makes it clear that the MS coupling is not expected to freeze 
within the 4-loop calculation for Nf < 5. 

In this situation we expect the BZ expansion to break down at Cq — 0.26 (cor- 
responding to Nf = 6). The coefficients of the BZ expansion are obtained directly 
from eq. (59), 

^ 11675 2 /145645559 5335 ^,^,\ 3 

x'fp = ao + ao^ + + C(3) oq^ (65) 

10272 V 17585664 428 V ^ ^ 

/ 92177206455497 587191201 \ 4 

^ V 1083839643648 ~ 13189248 ^ ^^'"V "'^ 

and numerically, 

xfp^ao + 1.1366 ao^ + 23.2656 ao^ + (-138.5630 + C4 -1) ao^. (66) 
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We see that the series is indeed ill behaved. For example, at oq = 0.26, the first 
three terms are roughly: 0.26, 0.0768, 0.409. In fact, the series seems to break down 
well above Nf ~ 6, but using only three terms it is difficult to estimate exactly 
where. 

We conclude that at least for Nj < 6 the MS coupling does not freeze, and 
that the "Landau pole" behavior of the 1-loop coupling in the infrared, persists in 
MS when higher-order corrections are taken into account. But is this conclusion 
true to all orders? We think that the answer is positive: since the BZ expansion 
in MS seems to break-down already at the order 0{al), it is hard to imagine how 
higher-order corrections could alter the situation. Still, it would be interesting to 
know how the BZ series behaves at higher orders. 

One approach [27] is to use the information that can be obtained on the MS 
f3 function at higher-orders is from the large- A^j limit. From the 0{1/Nf) terms 
calculated in [27] one can obtain all the C(l/ao) terms, denoted by (cf. (54) 
- (56)). In ref. [27] it was found that they are small, and that their effect on the 
value of the fixed point (and thus also on the question of its existence) is small. 

Another possibility is to use PA's [36] to estimate higher-order coefficients, and 
use them to calculate higher-order terms in the BZ expansion. Such an effort is now 
under way [37]. 

3.3 BZ Expansion for Time-like Quantities 

In Sec. 2 we found that the 'optimized scheme' approach does not give reliable results 
for the effective charge at freezing when applied directly to time-like quantities. As 
we show in Sec. 4, a consistent description of freezing for Re+e- or i?^ exists only 
when the space-like D-function freezes, and then the infrared limits are expected to 
be the same for all three quantities: 

L'(0)=i?e+e-(0)=i?r(0). (67) 

This is a direct consequence of the expected analyticity structure of the D-function, 
which indeed holds when perturbative freezing occurs. 

The 'optimized scheme' approach, in general, does not obey this requirement, 
since the terms that are related to the analytical continuation from space-like to 
time-like momentum, change the ECH (or PMS) (3 function (see eq. (28)). Fig. 
2 shows that the resulting difference between the values of the effective charges at 
freezing are insignificant above Nf = 8. However, for Nf ^ 8 these differences 
become larger. This can be interpreted as a sign that the three-loop results are not 
conclusively indicating freezing for these cases, which is true, but it is also related 
to the unjustified use of the ECH/PMS methods directly for time-like quantities. 

The BZ approach, on the other hand, yields the same expansion for the three 
effective charges, the D-function, Re+e- and Rr, just as one expects from (67). It is 
straightforward to check this: from the relations between the coefficients of Re+e- 
(27) (or R^- (31)) and those of the D-function, it is easy to obtain the relation between 
the corresponding coefficients of the different powers of Oq for the two quantities. 
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namely rij and dij (or Tij and dij). The substitution of rjj (or Tij) in terms of dij 
into the formulae for the BZ coefficients (eq. (63)), gives the same Wi coefficients 
as for the D-function. All the terms that result from the analytic continuation just 
cancel out between the different contributions, both for Re+e- and for Rr- 

3.4 BZ Expansion for various quantities 

In this section evaluate the coefficients in the BZ series for the fixed point in physical 
effective charges schemes for the various quantities discussed in Sec. 2.2. 

The calculation is straightforward using eq. (63) and the three loop coefficients 
from refs. [28, 29, 30, 33, 32]. Thus we go directly to the results, followed by a 
discussion. In the next section we study the implications of the Crewther relation 
between xo and xbj on the BZ expansion. Note that the 0{a^) coefficients cannot be 
explicitly calculated due to the lack of both the 5-loop coefficients of the P function 
and the 4-loop coefficients of the series for the different observables. 

1) For the polarized Bjorken sum rule: 

Bi 753 2 / 5930095 275 \ , 

= "° + 3424 "° + ll7585664 " 214 ^''^ 

/ 18602593666427 9470237 535 u \ ^ 

+ r 361279881216 ~ 137388" " "g" + "'-^ + ''''V "° 

where k^^ is the 4-loop Bjorken polarized sum rule coefficient at Nf = 16^, as 
in the general definition in eq. (61). The numerical results are: 

x% = ao + 0.2199 Cq^ - 1.2075 Oq^ + (-203.6939 + C4,-i + /ca.o) (69) 

2) For the non-polarized Bjorken sum rule: 

p^ _ 4589 2 /425842061 30809 805 \ 3 

"^^^ - ^° " 10272 + I 52756992 + "963~ ^'~l8 "° 

/ 38886699582523 327376259 2966545 \ 4 

^ V 361279881216 2198208 ^ 61632 ^ ^^'"^ ^ "^^'V 

or 

x^'p = ao - 0.4467 Oq^ + 0.1551 Oq^ + (-236.7461 + C4-1 + /s.o) Oo^ (71) 

3) For the GLS sum rule: 

rrs 753 o /239442925 3355 \ o 

a;^^* = ao H ai + (3 On (72) 

" 3424 ° V 52756992 321 ° ^ ' 

/_ 5706068695529 _ 535 171229447 ^ I \ a"^ 
^ V 120426627072 ~ ~8~ ~ 2198208 ^ ^"^'"^ ^ ^'V ^° 

or 

x^^-^ = ao + 0.2199 ao^ - 8.02495ao^ + (-210.3609 + C4,_i + ^3,0) Oo^ (73) 
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4) For the vacuum polarization D-function (and therefore, according to Sec. 3.3, 
also for Re+e- and Rr): 



^ _ 4177 2 /31250575 275 

X pp — Q/Q ~h — ~ — ~ do 



Cs) ao' (74) 



^FP . 3424 ' V 17585664 214 

/81595375713359 3893665183 2675 \ 4 

+ I 1083839643648 " 13189248 + ^ + "^'"^ + "° 
or 

= Oo + 1.2199 ao^ + 0.2323 oo^ + (-164.0075 + C4-1 + (^3,0) oo^ (75) 

5) For the static potential: 

y 8869 2 / 70824311 27 , 9 4 275 , \ o 

= ^° - 10272 1 17585664 + "64^ -214^^j"° ^^'^ 
/ 130549250005577 322462723 105075 4 

^ V 1083839643648 3297312 ^ ^^'"^ ^ ~ 219136 ^ 

315225 2\ 4 
+^7392"^ J 

or 

x^P = Oo - 0.8634 + 22.0945 + (-171.1353 + C4-1 + ^3,0) (77) 

6) For the derivative of the Higgs hadronic decay width: 

fr 31363 2 /486174653 275 \ o 

Xpp = ao H ai + C3 On (78) 

" 10272 ° V 52756992 214 ° ^ ^ 



/4675 , 294627948398435 982216871 , \ 

+ Cs Cs + C4 -1 o, 

V 48 3251518930944 9891936 ^ ' / 



xfp = ao + 3.0533 + 7.6707 + (-108.9778 + C4 -1) Oq (79) 
The conclusions are: 

a) The BZ coefficients for x^^, Xp^ and a;^^, are of order 1 up to 0{al). Esti- 

mation of the "radius of convergence" from the first few terms is difficult, but 
convergence of these series is not ruled out for any positive Nf. 

b) The BZ coefficients for XQ^^g, Xy^ and x^^ are relatively large. For XQ^g and 

Xy^ it is very difficult to estimate the "radius of convergence" from the first 
few terms since Wi coefficient is especially small, while W2 is rather large. For 
x^^ it seems that the series converges for ao ^ 0.36, corresponding to Nf ^ 2. 
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c) Clearly, for x^f^ is large due to the three loop light-by-light term that 
makes the only difference (cf. (20)) between the GLS sum rule and the Bjorken 
polarized sum rule, for which W2 is small. It may be surprising at first sight that 
light-by- light contribution which usually just a small correction to the three 
loop invariant p2 can make such a difference for the BZ expansion. The reason 
is basically the fact that the BZ expansion works around Nf — 16|, where 
the light-by-light term is not small, being proportional to Nf. In the vacuum 
polarization D-function, we neglected a similar light-by-light type term. As 
explained in Section 2.2, the correction expected in the D-function is smaller 
by a factor of {J2fQf)'^/Nf, compared to the GLS sum rule case. This term 
is therefore not expected to break the BZ expansion for the D-function. 

d) There is a noteworthy numerical cancelation between terms containing different 

irrational and transcendental numbers, that is responsible for the small W2 
values. This is particularly evident for the non-polarized Bjorken sum rule, 
where the rational term is roughly 8.07, the term proportional to (3 is 38.45, 
and the term proportional to C5 is —46.37, bringing the sum to about 0.155! 
Note that the case of the non-polarized Bjorken sum rule is special both in 
the fact that the 0{al) coefficient in the BZ expansion already contains a (5 
term, and in large cancelation between irrational numbers that occurs there. 

e) In x^^, which is a time-like quantity, all the tt^ terms that are related to the 

analytical continuation from space-like to time-like momentum that appear in 
the perturbative coefficients hi cancel out in the formula for the BZ series. 
This is in accordance with Sec. 3.3§. 



3.5 The BZ Expansion and the Crewther Relation 

In this section we study the consequences of the Crewther relation [21] (see section 
2.3) for the BZ expansions for the Bjorken sum rule and the vacuum polarization 
D-function. 

The first observation is that the Crewther relation can be used to evaluate the 
difference between the yet unknown 0{aQ) coefficients in the BZ expansion of the two 
observables. This observation is based on the fact that from the third relation in (37) 
one can obtain ^3^0 ~ c^3,0) the exact the difference between the 4-loop coefficients 
of the Bjorken sum rule and the vacuum polarization function at Nf = 16|: by 
substituting Po = the only unknown, namely S3, is eliminated. The result is: 

3439187 257719^ 535^ 
*,,„-rf3,c = ^^-^C3+ -a ^40.393 (80) 

Knowing /C3 — (^3,0 we look again at eqs. (69) and (75). If the BZ expansion 
for the location of the fixed point converges for both observables, then both O(ao) 

^Note that tt^ and tt^ terms do appear in the BZ expansion for Xy^ but these terms are not 
related to any analytical continuation. 
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coefficients sfiould be small. Thus also their difference should also be small. And 
indeed, using (80) this difference can be calculated directly, to give: 

- = _^^Z^ + ^ 0.70686. (81) 

^ ^ 35171328 107^^ ^ ^ 

Thus the Crewther relation is consistent with a very small O{ao'^) coefficient in the 
BZ expansions for both the D-function and the polarized Bjorken sum rule, with the 
caveat that it is also consistent with a common large contribution for both quantities. 

Another (related) observation is that thanks to the conformal Crewther relation 
at the fixed point (41), the BZ coefficients of the vacuum polarization D-function 
can be calculated directly from those of the Bjorken polarized sum rule, and vice 
versa. 

The BZ expansion of the D-function can be obtained from that of the Bjorken 
sum rule in two ways. In the first, one writes down the di coefficients in terms of 
the ki using (37). Then one uses eq. (63) to calculate the BZ coefficients, where one 
finds that all the Si terms just cancel out, to obtain: 

< = ^I' + lCf (82) 

The second way is by first calculating the BZ expansion for the Bjorken sum rule, 
using eq. (63), as in (68), and then substituting the entire series into the conformal 
Crewther relation (41), and finally expanding the rational polynomial again in terms 
of flo- The two methods are equivalent order by order. 

The practical implication of this is the following: if the BZ series indeed converges 
fast, then the rational polynomial one gets by substituting the BZ series for the D- 
function in the conformal Crewther relation (41) should itself be numerically close 
to the power expansion (74). Its higher-order coefficients should provide some rough 
estimate of the unknown higher-order terms. We present such an analysis in fig. 7, 
where we show both the straightforward BZ expansion for the vacuum polarization 
function, and the result of substituting the BZ expansion for the Bjorken sum rule in 
(41). The caveat is that the results of this analysis could be invalidated if there are 
higher-order effects that break the BZ expansion, and are common to the D-function 
and the Bjorken sum rule. 

3.6 The Derivative of the /3 function 

Following ref. [6] we now consider the derivative of the (3 function at the fixed point, 
defined by: 



dp{x) 

7 = 



dx 



-PqXfp 2 + Scxfp + 4:C2{xfp)^ + ■ ■ ■ (83) 



X=Xpp 
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Prom refs. [2, 6] we adopt the following form for the (5 function: 



^° + K + H{x) (84) 



(5{x) x"^ X ^{xpp — x) 

where H{x) is a power series in the coupling H{x) = Hq + HiX + ■ ■ ■, and ^ = j/Pq. 
7 is called a "critical exponent" since it determines the rate at which the coupling 
approaches the fixed point according to 

x-XFP^{QyAley (85) 

where A^q is the observable-dependent QCD scale [6] . 

It is well known that 7 is independent of the renormalization scheme, so long as 
the transformations relating the different schemes are non-singular (see ref. [35] and 
appendix B in ref. [6] and references therein). 

It is straightforward to obtain the BZ expansion for 7 by substituting eq. (58) 
in eq. (83) and expanding in powers of Oq. The result is: 

7 = oo (1 + giao + g2ao^ H ) (86) 

where 

91 = Ci,o 

92 = q - C3 _i - c|_i 

93 = cl o - 2C4 -1 - C3,o - 4C2 _i - 2C2,0C2 -1 

-4ci,oC3 _i - 5ci,oC2 _i - 6c2 -1C3 _i (87) 

The coefficients gi can be proven to be universal [15], in the sense that they are 
the same for any physical quantity, in agreement with the expectation that 7 is 
independent of the renormalization scheme in which the /3 function is defined. 
We now turn to the results. The BZ expansion for 7 gives: 

19 2 / 633325687 5335 \ o 

^ = "° + T"° + ll05513984 - W^^j "° + 
43834503808535 590624393, ^ \ 4 

+ ..n..o. C3-2c4,-0ao^ (88) 



270959910912 6594624 



and finally. 



7 = ao + 4.75ao^ - 8.89129ao^ + (269.43288 - 2c4,_i)ao^ (89) 



According to ref. [6], the divergence of the 7 series is not so bad to exclude the 
possibility of a perturbative fixed point for the physical case of = 0.36 {Nf = 2). 
We doubt this assertion, since the values the different terms in this case are: 1, 1.71, 
-1.15. Even for qq — 0.26 {Nf — 6) one obtains a rather slowly converging series, 
with the different terms contributing as follows: 1, 1.24,-0.6. 
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If we use Pade approximants (PA's) to estimate the unknown O(ao) term, we 
obtain a very inconsistent result. The [2/1] PA gives: 

gimp A ^ _ 192.494 (90) 

while the [1/2] PA gives 

g\V^]PA ^ 982. (91) 

The difference between the two PA predictions above provides an estimate of the 
uncertainty that we have assuming that g4 is close to zero. We therefore disagree 
with the assertion of [6] that a reasonable estimate for C4-1 can be obtained from 
the assumption that — 0. 

It is interesting that the gi-s are the same, not only for physical schemes, but 
also for MS , even though the MS coupling is unphysical. 



3.7 Conclusions 

The main question one would like to answer is at what Nf does the BZ expansion 
break down. Unfortunately, three terms in the expansion are not quite enough to 
provide a definite answer. In order to measure the convergence of the BZ expansion, 
we study the Nf dependence of the ratio of the C(ao'^) term and the partial-sum: 

W2ao^ /xpp provides some rough measure of the convergence of a series: a divergent 
series where all the terms are equal and positive yields a ratio of 1/3 above. If 
the signs oscillate, it yields 1. This ratio is presented in fig. 6 for the various BZ 
series. It is evident that the BZ expansion behaves differently for different physical 
quantities: while the expansion for value of the effective-charge at freezing seems 
to converge for any Nf for the polarized and non-polarized Bjorken sum rules and 
for the vacuum polarization D-function, it breaks down early for the GLS effective- 
charge (due to the light-by-light type terms), for the hadronic Higgs decay width 
effective-charge, for the static potential effective-charge and for the critical exponent 
7. It seems that the BZ expansion is reliable down to Nf — 12 in all cases, a point 
we shall come back to below. 

We conclude this section by comparing the picture one obtains for perturbative 
freezing from the two approaches studied, namely, finding zeros of the (3 function 
in an 'optimized scheme', and the BZ expansion. As a representative example, we 
choose the vacuum polarization D-function, and show in fig. 7 the value of the 
effective charge at freezing, as calculated by the ECH and PMS methods, together 
with the results from the BZ expansion. For the latter, we show both the result of a 
direct calculation (75) and the one obtained from the BZ expansion for the Bjorken 
sum rule using the conformal Crewther relation. 

We interpret the difference between the ECH and PMS results, as an intrinsic un- 
certainty of the 'optimized scheme' approach in this context. Similarly, we interpret 
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the deviation between the two BZ results as a measure of the intrinsic uncertainty 
of the BZ approach, related to the fact that we are using a power expansion in oq, 
rather that a more generic function of ao- 

From the comparison of the two approaches for calculating , i.e. ECH/PMS 
vs BZ, we conclude that for Nf 6, the three-loop result can lead to a perturbative 
fixed point - as shown in fig. 7, the two methods agree, and predict a relatively 
small effective coupling in the infrared limit: x ^ 0.3. We emphasize again that 
a zero in the truncated ECH/PMS f3 can be easily washed out by higher order 
corrections. An extreme example is provided by the Higgs decay width effective 
charge, for which at the 3- loop order it seems that there is a reliable fixed point, 
but in fact the perturbative series breaks down. 

The fact that various (space-like) effective charges run according to a very similar 
RG equation (at least up to 3-loops order) suggests that perturbative freezing will 
occur together and therefore that perturbative freezing at high enough order can be 
indicative for the existence of a fixed point in the full theory. 

The BZ expansions for the D-function, the polarized and the non-polarized 
Bjorken sum itiIcs show fast convergence up to order O(ao). The Crewther rela- 
tion is consistent with a small O(ao) coefficient for the D-function and the polarized 
Bjorken sum rule, but it is also consistent with a common large contribution at order 
C(ao) for both quantities. Prom the Crewther relation it is clear that if one of these 
quantities freezes, so does the other. 

An early break-down of the BZ expansion for physical quantities was identified 
for the critical exponent 7, for the static potential, for the derivative of the Higgs 
decay width and for the GLS sum rule. 

If we assume that existence of a genuine fixed point will be realized in a per- 
turbative manner, then we should expect convergence of the BZ expansion for any 
(infra-red finite) physical quantity. Using fig. 6, this leads to a prediction that 

Nf' > 12. (93) 

This result agrees with the results of Appelquist et al. [10], which are based on 
non-perturbative calculations and also with lattice simulations they refer to. On 
the other hand, it contradicts the results other lattice simulations [9]. 

Although this is outside the main subject of this paper, we emphasize again 
a nice feature of the perturbative expansions in QCD, that was noticed in two 
different occasions in the previous sections. This is the idea that the strong numerical 
cancelation between different irrational numbers (in QCD these are the Q terms) is 
usually not accidental and most likely provides an indication of some yet unknown 
deep relation that is encoded in the perturbative coefficients. Such a cancelation 
was found in the difference between the second RG invariants of the D-function and 
of the Bjorken polarized sum rule, as well as in the O{o^) term in the BZ expansion 
for the non-polarized Bjorken sum rule. 
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4 Analytic continuation and time-like quantities 



4.1 The D-function and Rp+ 



In this section we concentrate on the vacuum polarization D-function and the time- 
hke observables that are related to it through dispersion relations: Re+e- and Rr- 

The three loop analysis [1, 12], based on the (3 function in an 'optimized scheme', 
as briefly outlined in Sections 2.1 and 2.2, suggests that the D-function effective 
charge does not freeze for Nf < 4, while the related time-hke quantities do. The 
differences between the values of the time-like and space-like effective charges at 
freezing become significant already at Nf ^ 8. On the other hand, in the framework 
of the BZ expansion, the time-like and space- like quantities have the same expansion 
for the value of the effective charge at freezing. In the following sections we shall 
examine this issue on a deeper level. We will show that it is inconsistent to discuss 
perturbative freezing ofxji (or Xr) when the corresponding space-like effective charge 
has a "Landau-pole". We will also explain that the terms that are related to the 
analytical continuation are not supposed to change the value of the effective charge 
at freezing, i.e. xji{0) = Xt-{0) = xr){0), provided x oiQ) is well defined. 

We start by recalling [26] the analyticity properties of the D-function and the 
relations between the D-function and Re+e- ■ Ftom the optical theorem 

i?e+e-(s) = -Im{n(-s)} (94) 

TT 

where n(Q^) is defined in (21) and s > is a time-like momentum. From causality 
one expects that the only singularities of n((5^) are on time-hke axis, i.e. on the 
negative real axis: — —s with s > 0. The spectral density function, I3r{s) is 
defined by 

Us) - (95) 

Clearly Pr is also the P function for the coupling xr. Differentiating (94) one obtains 
a similar relation between Pr and the space-like D-function effective charge (22), 

f3Ris) = --lm{xDi-s)} (96) 

TT 

Based on the above, one can express the D-function as a dispersive integral over 

The relations between the coefficients of the corresponding effective charges xr and 
xd (27) arc directly obtainable from (97), as explained in ref. [26]. 
The inverse relation is: 

1 r+i^ , ,dii(-s') 1 r+i^ , ,D(-s') 

Re+e- is) = 7^ d'^TT^ = / ds'^^, 98 

2m Js-ie as' 2m Js-ie s' 
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where the integration contour hes in the region of analyticity of D{—s), that is, 
around the cut Re{s'} > 0. The contour can also be deformed to a circle, 

Relations (98) and (99) are also true for the corresponding effective charges, i.e. 
replacing Re+e- and D^Q"^) by xr and respectively, for instance, 

^r{s) -^.i ^xn{-s'). (100) 

2m J\s'\=s s' 

It is clear from eq. (100) that the infrared limit of Re+e-{s) equals that of the 
D-function if the latter exists: assuming the D-function does not have an essential 
singularity at s' = 0, we have xd{—s') — > Xd{0) and thus for small enough s the 
only singularity within the integration contour is the simple pole at the origin. Prom 
Cauchy's theorem one then obtains xii{s) — > xd{0)- 

From this argument one learns that the infrared limit of the exact Re+e- equals 
to that of the exact D-function. But does this hold in perturbation theory? The 
results of the 'optimized scheme' approach and the BZ expansion indicate that this 
is a delicate question. This issue is discussed further in the following sections. 

The analogous issue for R^ will not be discussed in detail in this paper, but most 
of our results are quite general and apply to it directly, since this quantity is related 
to the D-function in a similar way [31]: 

R^{ml) = — i — (l - 2a; + 2a;3 - x^) D{-xml) (101) 

Like in the Re+e- case, it can easily be shown that when D{—xml) — D{0) freezes, 
eq. (101) leads to an equality of the infrared limits Rt{0) — D{0). 

4.2 The analyticity structure of the D-function 

The physical quantity which can be measured directly is the time- like xr, or its 
derivative /3r(s), and not the space-like effective charge xdIQ"^). On the other 
hand, the perturbative calculation yields X£,{Q^). Naively, there is no problem to 
obtain (3^(3) from X£, via eq. (96). This is done by analytically continuing the 
perturbative xd{Q^) to time-like momentum = — s, where s > 0, and taking the 
imaginary part. Alternatively, one can obtain xr{s) from eq. (100). Unfortunately, 
the analyticity structure of perturbative xoiQ^) is inconsistent with this procedure. 
The relations between x^ and xr are based on the assumption that Xi:,{Q^) is 
analytic in the entire complex plane, excluding the negative real axis, < 0. 
On the other hand, a 1-loop perturbative result for xd{Q^) violates causality, since it 
has a "Landau-pole" at positive real Q'^ — A.qcd > ™ addition to the cut at < 0. 
In general, higher-order corrections to the jS function create a more complicated 
analyticity structure which also violates causality. Another way to see this problem 
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is that eqs. (98) and (99) lead to different results for xr with the same D function 
taken as input. This point is discussed in ref. [5]. 

In case of perturbative freezing, xoiQ"^) is finite for any positive real Q^. One 
might then be tempted to conclude that freezing saves causality. This is not nec- 
essarily the case, however, since in principle xd{Q'^) could still be singular at some 
complex Q^, while causality requires xd{Q^) to be analytical in the entire complex 

plane, apart from the time-like axis. Thus the rcsohition of the causality ques- 
tion requires full knowledge of the xj:,{Q'^) singularity structure. As we shall shortly 
demonstrate, the latter can be obtained explicitly. 

Before investigating in general whether freezing saves causahty, we study a spe- 
cific numerical example, where we demonstrate how perturbative freezing leads to 
an xoiQ"^) which is consistent with causality. We work with a 2- loop (5 function 

P{xd) = ^ = -PoXD^ (1 + c^d) (102) 

where t — ln((5^/Agjj), and for illustration purposes we take the hypothetical case 
I3q — 1 and c = —10. 

We will see that the solution of (102) defines a unique mapping from the entire 
complex Q2 plane (except the time-like axis) into a compact domain in the complex 
xd plane [38]. This domain does not contain the point xd = oo. 

A straightforward integration of (102) yields: 



/3oln(Q7Ay = - - — (103) 

X£) \ Xf) J 

In order to study the solutions of (103) in the complex plane, it is convenient to 
define: 

ln(QVA2fj)^p + ^r;, (104) 
where p > and — tt < < tt, and 

xd{Q'^) = r{]), Tj) + ik{p, Tj) = r + ik, (105) 

where both r and k are real. Eq. (103) can then be written as two equations for the 
real and imaginary parts 

1 , / 2 . 1 + 2cr 



and 



k r 1 

/^o^ = — o r7 + carctan k, c{r'^ + k'^) + r (107) 

+ L J 

where arctan[-y,, v] = arg{-u + iv} takes values in [— tt, tt). 

An infrared fixed point in the two loop /? function (c < 0) implies that for real 
(5^ > eq. (103) has a real solution < Xd < (— 1/c). In the notation of eqs. (106) 
and (107) this means that for 77 = and for any p > there is a solution with k — 
and < r < —1/c. Clearly, xd has a cut for real < 0, i.e. for r] — n. 
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In order to verify that there are no other singularities in the complex Q"^ plane, 
we explicitly find the domain in the complex xd plane into which the entire 
plane is mapped through (103). This is done by taking a contour around the cut 
and solving eqs. (106) and (107) numerically for A^g = 1. We choose the following 
contour in the plane: 

a) below the cut: = — s — ie, < s < oo, with e = 0.01; 

b) to the right of the cut: = Po + with po = lO"^*^ and — e < ^ < e; 

c) above the cut: = — s + ie with < s < oo. 

as shown schematically in the upper part of fig. 8. The resulting contour in the xd 
plane is presented in the lower part of fig. 8. Clearly, the particular choice of the 
contour around the cut is arbitrary, but a contour that is closer to the cut (smaller 
e and po) would correspond to an domain which is only slightly larger. We see 
that the entire complex plane is mapped into a compact domain m. Xd plane. 
Thus there are no spurious singularities in the complex plane and so causality is 
preserved. 

We emphasize that the solution of eqs. (106) and (107) described by this partic- 
ular mapping is not unique and there exist other branches of the solution. However, 
none of them corresponds to a real coupling xd along the positive real axis. It 
is this requirement that guarantees uniqueness of the solution. 

Returning to the more general case of a 2-loop /5 function, we will now show 
that perturbative freezing alone (c < 0) is insufficient to ensure that the analyticity 
structure of xd((5^) is consistent with causality. A further condition is required, 
namely that |c| > /Sq. 

The solution of the RG equation at the 2-loop order (102) in the complex 
plane can be written in terms of the so-called Lambert W function [39], which is 
defined by W{y) exp [W^(|/)] = y : 



xb{Q^) = -- 



z — — exp (—1 — Pot/c) = 

c ce 




(108) 



W{y) is a multi- valued function with an infinite number of branches, denoted by 
Wn{y)- We follow [39] as for the division of the branches and notation. The require- 
ment that xd{Q^) is real and positive for a real positive (at least for » -^eff)' 
is sufficient to determine the relevant branch: for c > the physical branch is 
W^i{y)^ taking real values in the range (— oo, —1), and for c < the physical branch 
is the principal branch, Wq{ii), taking real values in the range (—1, oo). 

We now need to check if the singularity structure oIxb{Q^) in (108) is consistent 
with causality. Let us first consider the case c > 0: W-i{y) has two branch points, at 
y = and at y — — 1/e, which are the endpoints of two cuts stretching to y = — oo. 
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The point at y = corresponds to 



oo. The point aiy — — 1/e corresponds 



to 



Qsing ~ -^eff" 



■c//3o 



^effexp(-c//?oln(c)) 



(109) 



which imphes that there is a "Landau singularity" on the space-hke axis. Thus for 
c > the singularity structure is inconsistent with causality. 

For c < 0, the physical branch W^iy) has only one branch point at |/ = — 1/e, with 
a cut stretching to y = — oo. In the complex plane this corresponds to branch 
points, i.e. "Landau singularities", given by (109). Two different possibilities exist: 

a) c < and |c| < id^. There is a pair of singularities in the complex plane. 

b) c < and |c| > /3o- There are no singularities in the first sheet and the 
pcrturbative solution is consistent with causality. In this case x^iQ'^^ maps 
the whole complex plane into a compact domain in the coupling plane. The 
latter situation was illustrated above in the example with /3o = 1, c = — 10. 

We will not consider here the higher-order /3 function. We expect that also in 
this case there are regions of parameter space resulting in an analyticity structure 
consistent with causality, i.e. allowing the mapping of the entire complex plane 
into a compact domain in the complex coupling plane. On the basis of the experience 
with the two-loop function, we conjecture that causality will be preserved by a 
wide subset, but not by the entire parameter space for which pcrturbative freezing 
holds. In general, the perturbative result does not have the analyticity structure 
implied by causality, and therefore does not obey the dispersion relation of (96) and 
(97). There is a way, however, to start from the perturbative /9 function and impose 
the dispersion relation on it: this is the so-called Analytic Perturbation Theory 
(APT) approach, which is the subject of the next section. 

4.3 Analytic perturbation theory approach 

The objective of the APT approach [5, 4] is to achieve the required analyticity 
structure, while retaining the correct perturbative behavior in the ultraviolet region. 

The technique is based on solving eqs. (106) and (107) on the negative real axis, 
i.e. for r\ — — tt. This can be done starting with the RG equation at any order. 
For instance, suppose we work at the 3-loop order and define the coupling a:((5^) in 
some arbitrary scheme and scale. Then xd = x ^ d\X^ -|- d2X^ . The imaginary part 
of xd{—s) yields the spectral density function (3r{s): 




TT 



(110) 



where x{—s) = r{t) +ik{t), and t = ln(s/Agg), where s is time-like, s > (c/. (105)). 
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The only case where the APT spectral function can be obtained in a closed form 
is for the 1-loop /3 function. Substituting c = and rj = — tt in eqs. (106) and (107), 
one obtains two algebraic equations for r(s) and k{s), whose solutions are: 

k{t) 



The spectral density at this order is simply 

Mt) - —kit) = - I (112) 

TT [(ln(s/A2)) +7r2 

Integrating the spectral density yields the time-like effective charge xr: 

XR{t) = arctan (n/t) + ^9 (-t) . (113) 

This is a positive, continuous and monotonically increasing function of s for any 
s > 0. Its infrared limit is l//9o- In the ultraviolet region this function approaches 
the 1-loop pcrturbative result xr{s) — > l/(/9ot). 

It has been recently emphasized (see the last ref. in [5]) that the results of 
the APT approach do not depend much on the renormalization scheme and scale. 
The reason for this stability (which has not been given in [5]) is rather simple to 
understand already from the 1-loop result: a scale or scheme transformation at 
this order amounts to a shift in A: A — > A. It is clear that the tt^ term in the 
denominator of (112) "hides" variations in A which are not too large. 

In order to study the higher order terms that make the APT result different from 
the standard perturbative 1-loop result, wc construct the expression for the spectral 
density (3^ (which is also the time-like f3 function) in terms of the effective-charge 
Xr. The perturbative regime corresponds to s > A^, so we can now ignore the 9 
function in the inversion of (113), obtaining 

t = 7^ ^ (114) 

tan{Po7^XR) 

We then substitute (114) in (112) and obtain: 

Po7r'\l + {tHPo7rxR)f 



Formally we work at the 1-loop level, so the higher-order terms in the r.h.s. of 
(115) are only a part of the full higher-order result. Still, we see that solving the RG 
equation on the time-like axis yields an infinite series of tt^/Sq terms. For instance, we 
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recognize the coefficient of tlie first correction, —P^n'^/S, as the difference between 
the ECH coefficients and pf" in eq. (28). Thus the APT approach can be viewed 
as a method of resumming the infinite series of terms associated with the analytical 
continuation. 

Following [5] we construct the corresponding space-like effective coupling, which 
is defined through the dispersion relation (97). The integral can be performed 
analytically, and yields 

^apt(<9^) = ^ 

The 1-loop APT effective coupling result (116) contains a first term that is just the 
1-loop perturbative result and a second term that exactly cancels the "Landau-pole" . 
Since this term is a power correction, it does not alter the perturbative ultraviolet 
behavior. 

Let us summarize the characteristics of the APT result, which are demonstrated 
above: by construction xjij)^{Q^) has a cut at Re{Q^} < and no other singularities 
in the complex plane, and in this sense it is appropriate to describe the D-function. 
Again, by construction, the APT result is consistent with the requirement that the 
time-like and space-like infrared limits are equal xji{0) — a;APT(0)- 1-loop the 
infrared limit value is xji^pj-{0) = l/f^o- 

We stress that there is an important difference between perturbative freezing, 
that leads to a finite infrared limit within perturbation theory, and the APT ap- 
proach, where xj^p^{Q'^) is not purely perturbative, as it contains power corrections 
that are due to the imposed analyticity. The non-perturbative nature of the APT 
result is less transparent on the time-like axis, where, as we saw, the APT approach 
can be viewed as a method to resum the infinite series of terms that are related to 
the analytical continuation. All-order resummations are in general dangerous. It 
is well known that the perturbative series itself has a zero radius of convergence, 
and that it is even non Borel-summable (see, for instance, [41]). Therefore, an all- 
order rcsummation of a partial series can yield any arbitrary result. One has to be 
convinced that a resummation procedure yields a result that is closer to the exact 
one, before utilizing it. It is important to understand that analyticity alone is not 
enough to set the infrared limit: one can still add further power corrections which 
can alter the infrared limit without violating the expected analyticity structure [40] . 

From the 1-loop APT result it is transparent that the terms in (3fi which are due 
to the analytical continuation have nothing to do with the existence of an infrared 
fixed point in QCD. Still, one can use the APT result to analyze the instabihty found 
in the predictions for the fixed point in the 'optimized-scheme' method (Sec. 2.5). 
As we saw, the fixed point of the "all-order" 1-loop APT result is at x^^ = I/Pq. 
On the other hand, the function /Sr of eq. (115), when truncated at order xj^, has a 
zero at x^^ = \/3/(7r/3o). However, when truncated at order x%, it does not have a 
non-trivial fixed point. At higher-orders, the j3 function has a non-trivial zero only 
when the series is truncated at an even power of x'^. Moreover, the convergence of 
the corresponding fixed point values to the limiting "all-order" value is quite slow. 



+ 



A^ 



ln(gVA2) A2-Q 



(116) 



34 



For instance, the deviation at order x% is 45%, at order x| it is 26%, and at order x)j 
it is 12%. This exercise shows that one cannot trust the value one gets by applying 
the 'optimized scheme' procedure to time-like quantities. 

Until now we discussed the APT approach for the 1-loop P function. It is in- 
teresting to see how the APT results change when higher loops are included. As in 
Sec. 4.2, we restrict our discussion to the 2-loop case. Starting from eq. (103) and 
performing analytical continuation by taking — —s with s > 0, we obtain: 

As mentioned in the previous section, inverting this implicit relation is non-trivial. 
In fact, there is an infinite number of solutions for xd{—s): 



c 1 + WniC^) 

(118) 

C± = --exp [-l-{t±t7r)Po/c] 
c 

for any integer n. As explained in detail in [42], uniqueness of the APT coupling is 
guaranteed only if the time-like solution X£){—s) is obtained from the physical space- 
like solution in such a way that for \Q'^\ ^ A^ the coupling xd{Q'^) is continuous 
for a continuous change of the phase of Q^. From here one proceeds, just as in 
the 1-loop case, by taking the imaginary part Pr = — (l/7r)Im{xz)(— s)} and finally 
calculating x^PT using (97), as in (116). 

While a full description of the analytic continuation from the space-like to the 
time-like axis using the various branches of the Lambert W function is postponed 
to [42], we briefiy present here some of the conclusions concerning the APT infrared 
limit. There are two possibihties: if c < — /Sq < (the perturbative solution is 
consistent with causality) the APT coupling coincides with the perturbative one 
and the infrared limit is a^APT(O) = xd{0) = — 1/c. Otherwise, for c > — /3o, 
^APt(O) = l//^o- Note the the latter is relevant both for —/Sq < c < cases, where 
there is a pair of complex singularities in the complex plane, and for c > cases, 
where there is one space-like singularity. 



5 Conclusions 

In this work we investigated the possibility that there is an infrared fixed point in 
QCD with Nf flavors, < A^j < 16|, that can be identifled from perturbative cal- 
culation. We examined the effective running coupling constant, defined from several 
QCD observables in the 'optimized scheme' approach and by the BZ expansion. 

We showed that several different (space-like) effective-charges behave similarly 
to one another. This suggests that freezing occurs for all of them together and 
therefore that perturbative freezing may be indicative of a genuine fixed point. 
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In general, the ECH/PMS approach, when apphed to space-hke quantities at the 
3- loop order indicates a possible perturbative fixed point for A^j ^ 5. It is clear, 
however, that the knowledge of higher-order corrections is essential for a conclusive 
answer. While for some observables the BZ expansion has small coefficients sug- 
gesting that it converges down to low Nf, in other physical schemes it breaks down 
quite early. Assuming that existence of a fixed point means that the BZ expansion 
should converge for any physical quantity, we get a prediction that Nj"^^* ^ 12. 

We emphasized a fundamental difference in QCD between the infrared behavior 
of quantities naturally defined for space-like and time-like momentum. We showed 
that perturbative freezing is a necessary (but not a sufficient) condition for consis- 
tency of perturbation theory with a causal analyticity structure. When the infrared 
finite perturbative coupling has a causal analyticity structure it coincides with the 
APT coupling. 

We emphasized that freezing leads to the equality of the time-like and space- 
like effective couplings at the fixed point. In this context it is important to stress 
that the terms due to the analytical continuation in the time-like (3 function are 
in principle not related to perturbative freezing. In the BZ expansion such terms 
cancel out, making it more reliable than the 'optimized scheme' approach, where 
they can lead to the appearance of spurious infrared fixed points. 
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Figure 1: The second renormalization- group invariant p2 as a function of the number 
of hght flavors for various space-hke (continuous hues) and time-hke (dashed hues) 
quantities. Also shown is C2, the second coefficient of the 13 function in MS (dotted 
line). Note the closeness of the curves for all the space- like quantities except the 
static potential {V). 
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Figure 2: The values of the effective charges at the fixed point, as a function of 
Nf, calculated in the ECH scheme (continuous lines) at the three loop order: upper 
line - the space-like vacuum polarization D-function, middle line - Rr and lower 
line - Re+e- ■ The dashed line represents the D-function effective charge at freezing 
calculated in the PMS scheme at the three loop order. Also presented is the ECH 
result for A^/ = 3 at the four loop order as calculated using predicted values for 
(see Sec. 2.5). 
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Figure 3: The differences between p2 values for various quantities and for the vacuum 
polarization D-function. Space-like quantities are plotted as continuous lines, while 
time-like quantities appear as dashed lines. The vertical axis is magnified 10 times 
compared to fig. 1. At small Nf all the space-like p2 invariants are close to each 
other. For Nf f» 16|, coincides with those of the related time-like quantities. 
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Figure 4: The renormalized difference TZ = J,- , „ as a function of for various 

values of Nf. The four lower lines correspond to Nf = (continuous line), Nf = 1 
(dashed line), Nf = 2 (dot-dash line), Nf = 3 (dotted line), and the four upper 
curves correspond to Nf = 4 (continuous line), Nf = 5 (dashed line), Nf = 6 
(dot-dash line), Nf = 7 (dotted line). 
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Figure 5: The renormalized difference TZ = L, /^J, as a function of Nf for various 

values of Nc. The four hnes correspond to iVc = 2 (continuous hne), Nc — 3 (dashed 
hne), Nf ^ A (dot-dash hne), Nf ^5 (dotted hne). 
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5 10 15 

Figure 6: The ratio of the 0{aQ^) term in the BZ expansion and the partial-sum 
{w2(io^ / Xpp) is shown for various quantities at the fixed point as a function of Nf. 
The dot-dashed fine represents the non-polarized Bjorken sum rule effective charge. 
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Figure 7: calculated from the ECH method (continuous line) the PMS 

method (dashed line) and the BZ expansion at the three loop order. For the BZ 
expansion, both the results of a direct calculation (dot-dash line) and a calculation 
that uses the Bjorken sum rule and the Crewther relation (41) (dashed line) are 
presented. 
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Figure 8: The mapping of the entire complex plane (except the negative real 
axis) into a compact domain in the complex coupling plane in the case of a 2-loop 
f3 function (102), with /?o = 1 and c = — 10. The upper plot shows schematically 
the contour around the cut in the plane as described in the text (Sec. 4.2). The 
lower plot shows the image of this contour in the complex coupling (xn) plane. The 
two fixed points (ultraviolet, at = 0, and infrared, at xd — — 1/c = 0.1) are 
denoted as squares. 
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